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Abstract 

The paper discusses extensions of the renormahzation group (RG) formalism for 3D 
incompressible Euler equations, which can be used for describing singularities develop- 
ing in finite (blowup) or infinite time from smooth initial conditions of finite energy. 
In this theory, time evolution is substituted by the equivalent evolution for renor- 
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Q , malized solutions governed by the RG equations. A fixed point attractor of the RG 

'c/3 I equations, if it exists, describes universal self-similar form of observable singularities. 

^^' This universality provides a constructive criterion for interpreting results of numerical 

O,. experiments. In this paper, renormalization schemes with multiple spatial scales are 
developed for the cases of power law and exponential scaling. The results are compared 

^^ ' with the numerical simulations of a singularity in incompressible Euler equations ob- 

K^ . tained by Hou and Li (2006) and Grafke et al. (2008). The comparison supports the 

^r I conjecture of a singularity developing exponentially in infinite time and described by 

^r^ ' a multiple-scale self-similar asymptotic solution predicted by the RG theory. 

^ ■ X 1 • 

O ■ 1 Introduction 

(N 

The question of whether the incompressible Euler equations in three dimensions can de- 
velop a finite time singularity (blowup) from smooth initial conditions of finite energy is the 
rS ■ long-standing open problem in fluid dynamics [11 [21 [3l H] . This question is of fundamental 

c^ . importance, as the blowup may be related to the onset of turbulence [5] and to the energy 

transfer to small scales [6]. Direct numerical simulations provide a powerful tool to probe 
the blowup. However, despite of large effort, we are still far from having a definite answer. 
There is a number of blowup and no-blowup criteria, which are useful in numerical 
simulations to detect a finite-time singularity. The widely used criterion is due to the Beale- 
Kato-Majda theorem [7], which states that the time integral of maximum vorticity must 
explode at a singular point. Several criteria, which also use the direction of vorticity, are 
developed by Constantin et al. [8], Deng et al. [HllIO] and Chae [H]. See also [H [121 [13] ^^^ 
references therein. 
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The history of numerical studies is summarized in [3l HSj [Tl]. Conclusions based on 
numerical simulations vary depending on initial conditions and numerical method. However, 
none of the results seem to be sufficiently convincing so far. Apparently, the success of 
numerical simulations requires further development of the theory. 

In this paper, we consider the renormalization group (RG) approach to the blowup prob- 
lem. The RG method is famous to capture sophisticated critical phenomena characterized 
by scaling universality, e.g., critical phenomena in second-order phase transitions [T5] and 
period-doubling route to chaos [12], see also the review in [T7]. There are various applica- 
tions of this method to problems of fluid dynamics [18] . Universality of self-similar blowup 
was explained in [19] for inviscid shell models of turbulence using the method, which is sim- 
ilar in spirit to the RG approach. Analogous universal self-similar blowup was observed in 
cascade models of the Euler equations [2D]- The RG approach to the blowup problem for 
incompressible Euler equations was discussed in [2T1 [221 123] in the case of a single spatial 
scale. In this paper, we study extension of the RG method for multiple spatial scales in the 
cases of power law and exponential scaling. Note that the universality predicted by the RG 
theory provides a constructive criterion for interpretation of numerical results. 

We start by illustrating an idea of the RG method on the inviscid Burgers equation, 
where the blowup phenomenon is simple and well known. Here the RG equation is derived 
for solutions renormalized near a singularity. Fixed points of the RG equation correspond 
to self-similar blowup solutions. Existence of an attracting flxed point explains universal 
scaling of the blowup [TH |23] . 

The RG equations for incompressible Euler equations are derived flrst for the case of 
a single spatial scale. Fixed points of the RG equations correspond to exact self-similar 
solutions. A part of this theory corresponding to flnite time singularities with power law 
scaling was considered in [2X11221 [23] . We extend the RG formalism to the case of singularities 
developing exponentially in inflnite time. Note that exponential scaling in this problem was 
suggested in [2U [2S] based on numerical results. 

The main contribution of this paper is the development of the RG formalism with multiple 
spatial scales. Self-similar solutions with such scaling cannot be exact solutions of the Euler 
equations. However, they may serve as asymptotic solutions. This is proved by introducing 
a special renormalization of the pressure term in the RG equations. Attracting flxed point 
solutions of the RG equations, if they exist, describe universal asymptotic form of observable 
flow singularities. Two types of scaling are considered, which correspond to flnite time 
(blowup) and inflnite time (exponential) singularities. 

The asymptotic forms of singularities provided by the RG theory are tested using the 
results of numerical simulations obtained by Hou and Li [211 [2Z] and Grafke et al. [21] . The 
comparison supports the conjecture of [2H [22] on exponential scaling of a singularity, which 
is consistent with a multiple-scale self-similar form of solution allowed by the RG theory. 

The paper is organized as follows. Section [21 describes the RG theory for the inviscid 
Burgers equation. Section [3] describes the single-scale version of the RG formalism for the 
incompressible Euler equations. Section [H develops the RG theory in the case of multiple 
scales. Section [5] extends the results to the case of multiple-scale exponential singularities. 
Section [B] compares the theory with known numerical results. Conclusion summarizes the 
contribution. 



2 Attractor of renormalized Burgers equation 

In this section we demonstrate the idea of the RG approach on a simple example of the 
inviscid Burgers equation 

ut + uu^ = 0, (2.1) 

which has the well-known classical solution leading to a singularity in finite time (blowup). 
This example contains many features of the RG formalism for the incompressible Euler 
equations developed in the next sections. 

Solution of equation 02. ip is given implicitly by the method of characteristics as 

u{x,t) = uo{z), x = z + Uo{z)t, (2.2) 

where uo{x) is a smooth initial condition at t = and z is an auxiliary variable. The blowup, 
u^ — 7- oo, occurs when dx/dz = 0. Expressing this derivative from the second expression in 
(12. 2p yields 1 + dzUo{z) t = 0. Therefore, one finds the time tf, and position Xb = z oi the 
blowup from the condition 

U = min (-Tr\^) ■ (2-3) 

We assume that ut, = u{xh,tb) = at the blowup point. This condition can always be 
satisfied by means of the Galilean transformation, which is a symmetry of (12. ip . 
Following [TS], we introduce the variables 

x' = X - Xb, t' = tb- t, (2.4) 

where t' > is the time interval to the blowup, and consider the renormalized solutions 
defined as 

uix,t)=&~'^U{^,T), i = x'/t'\ T = -\ogt'. (2.5) 

In the new variables, the blowup corresponds to r — )■ oo. Equation for f/(^, r) is found from 
dnD, (EH) as 

U, = [a- l)U - aiUf. - UU^. (2.6) 

For a = 3/2, this equation has a stable fixed point solution U{^, r) = t/*(^) given implicitly 
by US] 

e = -t/* - CUl C > 0. (2.7) 

Moreover, one can show [22] that 

limt/(e,r) = [/,(0 (2.8) 

for any blowup solution with nondegenerate minimum (12. 3p and some constant C. Note that 
this result extends to general scalar conservation laws [23] . 

Expressions (12.50 with a = 3/2 and (12.80 imply the asymptotic relation 



u 



{x,t)^t'^/^U.{x'/t'^/^), x'^t'^'\ t'^0. (2.9) 



An important observation is that this asymptotic relation is valid in the vanishingly small 
neighborhood x' ~ t'^/"^ of a singularity. For large ^ = x'/t'^^^, using the approximation 
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Figure 1: (a) Solution of the inviscid Burgers equation at t' = 1, 1/2, 1/4, and 1/8, where 
t' = tf, — t is the time to blowup. Renormalization boxes, which scale as x ~ t'^'"^ and 
u ~ t'^/^, are shown, (b) Renormalized solutions converge to the universal function U^{^), 
which is shown by the bold red curve. 



U^: ~ —{^/Cy^^ following from (12.71) in (12. 9p leads to the classical cubic- root expression 
u{x,t) ^ —(x/Cy^'^. Therefore, the wave profile u{x,t) contains a universal self-similar 
"core" of the form (12. 9p shrinking to a point as the blowup is approached. This "core" has 
the cubic root x-dependence for x' ^ t'^^'^, where it is "glued" to the rest of the solution. 
An example of convergence to the fixed point solution U^ is shown in Fig. [H 
The scaling symmetry 

u{x,t) \-^ bou{x/bo,t) (2.10) 

can be used to set the value C = 1 in (12. 7p . The choice of tb, Xb at the blowup and the 
condition Ub = are important for the convergence to f/* in (12.80 . as small errors in satisfying 
these conditions lead to instabihty [18]. The values of Xb, h and Ub can always be adjusted 
by using the symmetry transformations 

u{x,t) \-^ u{x + bi,t), u{x,t) \-> u{x,t + b2), u{x,t) \-^ u(x + b3t,t) — bs. (2.11) 

In summary, the blowup in the inviscid Burgers equation is related to the evolution of the 
renormalized function f/(^, r) near the fixed point 6^*(0 in functional space demonstrated 
schematically in Fig. [2l The fixed point f/* has a stable manifold Ms of codimension 4. 
The four extra dimensions correspond to solutions related by symmetries (12.101) and (12. lip , 
which determine 4-dimensional surfaces E of equivalent solutions. Surfaces E intersect Ms, 
so that one can move any point to the stable manifold by the symmetry transformation. 
This, in turn, implies universality of the limiting renormalized solution. 

One can check that the shift r i— )■ r -(- At in (12. 5p induces renormalization in space-time. 
This shift can be viewed as an action of the renormalization operator 



TZat 



f/(e,r)^^f/(e,r + Ar). 



(2.12) 



Since the right-hand side of (12.61) does not depend on r explicitly, this operator generates a 
one-parameter renormalization group with 7lri+T2 = T^nT^Ti- Existence of such a group is 
an important property underlying the above construction, similarly to other applications of 
the RG theory [El [16l [I?] . 




Figure 2: Schematic local structure of the fixed point U^, of the RG equation (12.61) . Ms 
represents the stable manifold. Points of each surface E correspond to solutions related by 
symmetries fl2.1Up and fl2.1ip . Universal structure of the blowup is explained by the RG 
dynamics on the stable manifold. 

3 Exponential renormalization of incompressible Euler 
equations 

The Euler equations governing the fiow of ideal incompressible fiuid of unit density in three 
dimensional space x = {xi,X2,X3) are 

Ui + u ■ Vu = -Vp, V ■ u = 0. (3.1) 

We consider solutions u(x, t) with smooth initial conditions of finite energy, J d^xWuW^ < oo. 
First, let us assume the finite time blowup, when the fiow forms a singularity at (xb,^;,). 
Using the Galilean transformation, one can set u;, = u(xb, U) = 0. Leray [30] suggested to 
consider self-similar solutions of the form 

u = t'("-i)U,(x7t"'), p = t'(''^-2)p,(x7t"*), a>0, (3.2) 

in the study of blowup, where 

x' = X — Xfe, t' = tf, — t. (3.3) 

Finite energy solutions of this form cannot be realized globally [3T| [1] . However, the global 
existence is not required in the RG theory, since the self-similar expression is valid asymp- 
totically in a vanishingly small neighborhood of a singularity, see (12. 9p . 
The renormalized solution is defined similarly to (12. 5 p as 

u(x,t)=t'('^-^)U(^,r), p(x,t)=t'(2'^-2)p(^,r), ^ = x7t'^ r = -logt'. (3.4) 

Substituting these expressions into (13. ip . one obtains the RG equations (see [211 [22]) 

Vr = {a- 1)U - (a^ + U) ■ V^U - V^P, V^ • U = 0. (3.5) 

A fixed point solution U*(^), -P*(^) satisfies the equations 

(a - 1)U, - {a$, + U*) ■ V^U* - V^P* = 0, V^ ■ U* = 0, (3.6) 



Substituting these expressions into fl3.ip . one obtains the RG equatic 



and determines the self-similar solution (13. 2p . The vorticity computed for the velocity (13. 2p 
grows as a; ~ 1/t', which agrees with the Beale-Kato-Majda theorem [7]. 

Fixed point solutions describing the blowup are not known, but something can be said 
about their stability assuming that they exist [22]. Stable fixed points would describe ob- 
servable blowup phenomena. Different types of attractors of the RG equations like, e.g., 
periodic solutions are also relevant for the blowup problem [32] . 

The RG formalism can be extended to the case of a singularity developing exponentially 
in infinite time. Assuming that the solution is regular for all t > (no finite time blowup), 
we consider the renormalized functions U(^, t) and -P(C, t) defined as 

u(x,t) = e-^*U(^,t), p(x,t) = e"2^*P(^,t), | = x7e-^ 6 > 0. (3.7) 

tions 

Ui = 6U - (6| + U) ■ VgU - V^P, V^ ■ U = 0. (3.8) 

A fixed point solution U*(£), P*(^) satisfies the equations 

6U* - {hi + U,) ■ V^U, - V^P* = 0, Vj ■ U, = 0, (3.9) 

and determines the self-similar solution 

u(x, t) = e-'*U(x7e-'*), p(x, t) = e-^'^Pi^ /e-"") (3.10) 

describing a singularity developing exponentially as t -> cx). Note that one can always set 
6 = 1 by time scaling. 

The vorticity computed for the velocity (I3.10p remains finite, but second spatial deriva- 
tives of the velocity grow exponentially with time. Since vorticity grows in all numerical 
simulations of the incompressible Euler equations, it is unlikely that stable fixed points 
of (13. 8p exist. However, analogous exponential singularities with multiple scales described 



below seem to be good candidates for describing observable phenomena. 

We remark that equations (I3.7p - (l3.10p follow from (I3.2p - (l3.6p in the limit a — ^ oo. For 
example, the fixed point equations (13. 9p are obtained from (13. 6p in the limit a — )■ oo after 
multiplication by fe^/a^ and the substitution U* = (a/6)U*, P^, = (a/6)^P*. Similar relation 
is established for other expressions by taking t' = 1 — bt/a so that 

t'^ = (1 - bt/a)'' -^ e-^ r = - logt' ^ bt/a (3.11) 

as a — 7- oo. 

Recall that, in the incompressible Euler equations, the pressure is determined by the 
velocity field at the same time. The same, of course, is valid for the RG equations. Indeed, 
computing divergence of both sides of the first expression in (13.80 . the left-hand side vanishes 
due to the incompressibility condition and the right-hand side yields Poisson's equation for P. 
Its solution is well determined if the vector field U(^, t) decays sufficiently fast as ||^|| — )■ oo. 
However, following the results of Section |51 we can expect that the fixed point solution 
U*(^) of the RG equations is unbounded for large ||^||. In this case the problem (13. 9 p for a 
fixed point is not well-posed. This inconsistency is removed in the multiple-scale RG theory 
considered in the next sections. 
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4 Multiple-scale RG formalism 

Self-similar expressions of the form (13. 2 p or f l3.10p describe singularities with a single spatial 
scale, a; ~ t'" or X ~ e~^*. On the other hand, numerical simulations of incompressible 
Euler equations available in the literature demonstrate very thin singular structures implying 
existence of at least two different scaling laws. For example, the two scales proposed in [33] 
are xi ~ 1/Vt^ and X2 ~ 1/t'. In this section, we generalize the RG formalism for the case 
of multiple-scale self-similar solutions. 

Let us assume that solution u(x, t) with smooth initial condition of finite energy blows 
up in finite time at (xb,tb). Using the Galilean transformation, we set Uf, = u(xfe,tfc) = 0. 
Let us introduce the diagonal matrices 

t'^ = diag(t"^^t'»^r^), A = diag(ai,a2,a3), (4.1) 

which generalize the scaling x' ~ t'"- of Section [3] to the multiple-scale case. In what follows, 
we assume 

< ai < as < as (4.2) 

with a single dominant power a^. The case of equal 02 and 03 can be analyzed similarly. 
The renormalized solution is defined as 

u(x,t)=t'(^-i)U(|,r), p(x,t)=t'(2'^^-2)P(^,r), | = t'^^x', r = -logt', (4.3) 

where I is the identity matrix. The first relation written for each vector component has the 
form 

«^. =t'K-i)f/^.(^,r), J = 1,2,3, ^ = (xl/t'»^X2/t"^^X3/t'"^), (4.4) 

where three different scales are used for different space directions. 
Substituting (14.31) into (13.11) yields 

U, = (A - I)U - (A^ + U) ■ V^U - CV5P, V5 • U = 0, (4.5) 

where C is the diagonal matrix 

Taking time derivative of (14.61) yields the equation 

Cr = 2(A - a3l)C. (4.7) 

with the initial condition 

C(0) = I. (4.8) 

We see that equations (14. 5p . (14. 7p do not depend explicitly on the renormalized time r, 
which is the key point of the RG approach. Thus, (14. 5p and ( 14. 7p determine a renormalization 
group parametrized by r, as it is explained in the end of Section O 

Let us analyze fixed points of the RG equations. According to (14.61) and (14. 2p . we have 
the constant solution 

a = diag(0,0,l), (4.9) 
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which is a fixed point attractor of (14.71) with initial condition (14. Sp . Then, the fixed point 
solution U*(^), P*{$,) of (14.51) is determined by the equations 

(A - I)U, - (A^ + U,) ■ VgU, - aV^P. = 0, V^-U, = 0, (4.10) 

According to (14. 3p . the fixed point solution defines the self-similar fiow 

U(x,t) = t'(^-^)U,(t'-^x'), p(x,t) = t'(2a3-2)p^(^/-A^/)^ ^4^^^) 

For each velocity component, the first expression reads 

w,(x,t) = t'^'^^-'^U,,ix[/t"^\x'Jt'^^xyt'^^), J = 1,2,3. (4.12) 

Velocity field and pressure in (14. lip are not exact solutions of the Euler equations (13. ip . 
since multiple scaling is not a symmetry of the inviscid incompressible fiow. However, for 
solutions U attracted to the fixed point U*, (14. lip yields the asymptotic form of the fiow. 
This asymptotic expression is valid in the vanishingly small neighborhood of the blowup, 
i.e., for 

x[ ~ t"'\ x'^ ~ r^ x'^ ~ r^ t' -^ 0, (4.13) 

corresponding to constant values of ^. This fact can also be checked directly by substituting 
(14. lip into (13. ip and noting that the pressure terms of the first two Euler equations are 
asymptotically small. The asymptotic blowup solution (14. lip is observable, if (U*,P*, C*) is 
an attractor of the RG equations (allowing for irrelevant unstable modes related to system 
symmetries, e.g., space translations and rotations). 

The vorticity computed for the velocity (I4.12p grows as oj ~ f'iai-a^-i) (determined by 
the dominant derivative dui/dx^. Since a^ > Oi, this agrees with the Beale-Kato-Majda 
theorem [7]. Note that the first two components of the vector field (I4.12p satisfy exactly the 
Euler equations with vanishing pressure, i.e., we have 

dui/dt + u-Vui = 0, du2/dt + u-Vu2 = 0, V ■ u = 0, (4.14) 

as one can check by the substitution of (I4.12p into (I4.14p and using (14. 9p . (I4.10p . The Euler 
equation for the third component remains unchanged and determines the pressure. We see 
that the pressure decouples from system (I4.14p for the velocity components. Hence, the 
inconsistency related to determining the pressure function mentioned in the previous section 
does not appear in the multiple-scale theory. 

5 Multiple-scale exponential singularity 

In this section we consider extension of the multiple-scale RG theory to the case of exponen- 
tial scaling introduced in Section [HJ Let us assume that a fiow u(x, t) with smooth initial 
condition of finite energy is regular for all t > and forms a singularity as t — > oo at the 
point Xfe with u^ = 0. Exponential renormalization with multiple scales is defined using the 
diagonal matrices 

e^* = diag(e^^*,e*^*,e^3*), B = diag (6i, 62, ^s)- (5.1) 



We consider the case 

0<bi<b2<h (5.2) 

with a single dominant exponent 63. The case of equal 62 and 63 can be analyzed similarly, 
and the single-scale solutions of Section [3] correspond to fci = 62 = &3- 
The renormalized solution is defined as 

u(x,t) = e-«*U(^,t), p(x,t) = e-2''3*P(^,t), ^ = e^V. (5.3) 

The first relation written for each vector component has the form 



u 



'] 



e-'^%i^,t), J = 1,2,3, ^ = (xi/e-^i*,X2/e-^^*,X3/e-^^*). (5.4) 



Substituting fl5.3p into (13.1 p yields the RG equations 

Ui = BU - (B| + U) • VgU - CV^P, Vg ■ U = 0, (5.5) 

where C is the diagonal matrix 

C(t) = e-^'^'e^^' = diag ^e^t(b,-b^) ^^2tib,-b,)^ ^y (5 g^ 

This matrix satisfies 

Ct = 2(B - 63l)C, C(0) = I. (5.7) 

As in the previous section, we obtained the consistent RG theory, where equations (15. 5p . 
(15. 7p do not depend explicitly on time. Equations (15.70 have the attracting constant solution 
(14. 9p . Thus, the fixed point solution U^=(^), P*{$,) of (15. 5p is determined by the equations 

BU, - (B^ + U,) ■ V^U* - CVgP, = 0, Vg ■ U* = 0. (5.8) 

According to (15. 3p . the fixed point defines the self-similar asymptotic solution 

u(x, t) = e-^*U,(e^V), p(x, t) = e-^'^^P^e^'x'). (5.9) 

For each velocity component, the first expression reads 

Mj(x,t) = e-^^'U,j{x[/e-^'\xye-''^\xye-^''), 3 = 1,2,3. (5.10) 

The vorticity computed for the velocity (I5.10p grows as cj ~ eibs-bi)t cletermined by the 
dominant derivative dui/dx^. 

Similar to (14. lip , velocity field and pressure in (15. 9p are not exact solutions of the Euler 
equations (13. ip . but they provide an asymptotic solution valid in the vanishingly small 
neighborhood of a singularity 

x\ ~ e-^^\ 4 ~ e-^^\ X3 ~ e"''^*, t -^ cx). (5.11) 

This solution describes an observable singularity, if the fixed point (U*, P*, C*) is an attractor 
of the RG equations (allowing for irrelevant unstable modes related to system symmetries). 
In this case, (15.100 is a universal asymptotic form of a singularity. Note that the coefficients 
bi, 62 and 63 are multiplied by the same positive factor under a change of time scale and, 
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Figure 3: Dependence of the maximum vorticity in log-scale on time reconstructed from 
(a) Fig. 9 in [27] and (b) Fig. 1 in [2B]- One can notice the asymptotic exponential time 
dependence (straight dashed lines). 

hence, only the ratios 61/63 and 62/63 are expected to be universal. As in Section [3l one can 
obtain the multiple-scale exponential singularity expressions as the limit of the formulas for 
blowup in Section m with Qj = bja and a — )■ 00. Also, the asymptotic solution (15. 9 p is the 
exact solution of (I4.14p . as it can be checked by the substitution using (15. 8p . 

We remark that a scaling of pressure different from (14. 3 p and (15. 3p can also be considered 
in the RG scheme. For example, the scaling p(x, t) = P{$,, t) with 61 = 62 = 1 was suggested 
in [25]. For such scaling, however, the problem for fixed points of the RG equations is not 
well-posed, due to similar reasons as described in Section El 



6 Self-similar exponential singularity in numerical sim- 
ulations 

In several studies, numerical solutions were interpreted in favor of finite-time blowup, sug- 
gesting the growth of maximum vorticity as max|a;| ~ 1/t', see, e.g., [311 ESj- However, 
this asymptotic relation was not confirmed in later computations [36l[2H|- Several numerical 
studies [251 EH ISSj suggested exponential time dependence of vorticity. Self-similarity of 
numerical solutions was discussed in [2^ l25l 133). See also the review of numerical results in 



Fig. El shows dependence of maximum vorticity on time reconstructed from Fig. 9 in 
[27] and Fig. 1 in [28]. Initial conditions in these simulations have the form of antiparallel 
vortex tubes in [27] and of Kide-Pelz 12 vortices in [28]. Logarithmic scale is used for the 
maximum vorticity. One can see that max |cj| is approximated very well by the exponential 
time dependence (straight dashed lines) for large times. Similar exponential behavior was 
observed along a Lagrangian trajectory passing near the singularity, see Fig. 5 in [25] . 

In this section, we use the RG formalism to test the conjecture [25] on exponential scaling 
and self-similarity of flow near a singularity. The appropriate self-similar solution is described 
by (I5.10p . For maximum vorticity, it yields the asymptotic relation max|cij| ^ uqc 



{b3-bi)t 
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given by the dominant term dui/dxs, in agreement with Fig. [3l 
It is convenient to use the Fourier transformed form of fIS.lOp as 

uj{k,t) = e-^'^+''+'^+'-'^%,{h/e''^\h/e'^\h/e''^'), j = 1,2,3, (6.1) 

where u(k, t) is the Fourier transform of u(xb + x', t) and U*(k) is the Fourier transform of 
U*(^). The asymptotic relation fIS.lOp is vahd in a small neighborhood of a singularity in 
physical space (15. lip . Hence, the self-similarity must be observed for large kj, namely, for 

kj ~ e^^\ J = 1,2,3, t^oo. (6.2) 

In particular, it follows from (16.21) with the conditions (15. 2p that, for large t, 

ki < k3, k2 < /C3, k= ||k|| ^ ks. (6.3) 

This means that the self-similar part of solution is concentrated near the axis k-s in the 
Fourier space. Expressions (I6.1l) - (l6.3p imply the self-similar asymptotic behavior of the 
energy spectrum 

E{k, t) ^ e-^^'E{k/e^''), k ^ ka ^ e^'\ t -^ cx), (6.4) 

which scales as the dominant variable ks. The function E{k) and the constant bs can be 
expressed in terms of U^,(k) and bj. 

Below we verify the self-similarity hypothesis by checking (16. 4p for the numerical data 
reconstructed from Figs. 9 and 10 in [26]. Fig. |4](a) shows the graphs of energy spectra 
E{k, t) at times t = 15, 16, 17 before and after the scaling chosen as 

E{k,15), e^-^E{e'^-^'^k,16), e^Eie'^-^^k,!?). (6.5) 

Good matching of the scaled profiles is observed. The scaling exponents change almost 
linearly in time in agreement with the asymptotic expression (16.41) . Similar behavior was 
observed in p^ [25] based on the approximation E{k,t) = c(t)A;~"'^*)e~^'^*^*^'^. 

More detailed analysis can be carried out using the spectra Qn.{k,t) = k'^^E{k,t) of 
generalized enstrophies considered earlier in [21]. According to (16.41) . these spectra must 
scale exponentially in time as 

nn{k, t) ^ e^^'''"-'"'^'nn{k/e^''), k ~ e'"', t -^ oo, (6.6) 

where ^n{k) = k'^"'E[k). The reason to consider Qn{k,t) instead of E{k,t) is the following. 
For sufficiently large n, the factor A;^" (corresponding to the spatial mean-square derivative) 
suppresses the energy spectrum in the region of small wave numbers, k <C e^^*, where the 
asymptotic relation (16.61) is not valid. As a result, the function VLn{k,t) is large for k ~ e^^^, 
where the self-similarity is expected, and gets small far from this region. 

Fig. m^b) shows the function VtQ{k,t) divided by its maximum value at times t = 8, 10, 
12, 14, 15, 16, 17. The shape of this function is almost independent of time, as shown 
in Fig. m^c), where the proffies are shifted to match at the maximum. Fig. [3] presents 
the values of max fig and the corresponding values of k for the proffies in Fig. ID^b). The 
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Figure 4: Energy spectra reconstructed from Figs. 9 and 10 in [26]. (a) Energy spectra before 
and after scaling (16. 5p . Log-scale is used for E. The scaled spectra match in agreement with 
(16. 4p . (b) Normalized spectra QQ{k,t) at different times with log-scale for k demonstrating a 
"traveling wave" moving toward large k. (c) Profiles of the previous figure shifted to match 
at the maximum, demonstrating self-similar asymptotic form in agreement with (16. 6p . 

plots demonstrate asymptotic exponential dependence for t > 14 (straight dashed lines 
in logarithmic scale). Fig. [5] can be compared with the asymptotic exponential growth of 
maximum vorticity in Fig. [3t^a), corresponding to the same numerical simulation. Note 
that the blowup in the inviscid Burgers equation (12. ip as well as in inviscid shell models of 
turbulence is characterized in the Fourier space by the behavior very similar to Fig. Hl^b) 
under appropriate renormalization, see [T9| [23]. 

The presented analysis of numerical results supports the conjecture that the inviscid 
incompressible flow has a singularity developing exponentially in time and having asymptotic 
self-similar structure (I5.10p . The RG theory of Section [5] suggests that the renormalized 
velocity field and ratios of scaling coefficients may be universal, i.e., independent of initial 
conditions. This universality is a powerful criterion, which can be checked numerically. 
Such a test, however, is nontrivial, since the universality may be sensitive to symmetry 
transformations, see Section [2l 

7 Conclusion 

The problem of existence and structure of singularities developing in finite time (blowup) or 
infinite time from smooth initial conditions of finite energy in incompressible Euler equations 
is considered. These singularities may be studied using the renormalization group (RG) 
approach. The central point of this approach is deriving the RG equations, which determine 
the fiow evolution combined with renormalization of space, time and velocities. A fixed 
point attractor of the RG equations, if it exists, describes a universal self-similar form of 
observable singularities. 

In this paper, we described extensions of the RG formalism, which include multiple scales 
and different scaling laws. We explained the possibility of two types of universal self-similar 
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Figure 5: Analysis of energy spectra reconstructed from Figs. 9 and 10 in [26]. Shown 
are the time dependence of (a) m.a.xQQ{k,t) and (b) the corresponding values of k for the 

k 

profiles in Fig. Hl^b). Log-scale is used for the vertical coordinate, demonstrating asymptotic 
exponential dependence (dashed lines). 



flow structures valid asymptotically in a small neighborhood of a singularity. The first type 
describes formation of a finite time singularity (blowup) with the power law scaling x' ~ t"^^ , 
j = 1,2,3. In the limit aj — )■ oo, we obtain the second type corresponding to solutions 
with exponential scaling, x' ~ e"^^*, which describe singularities developing exponentially 
in infinite time. Such a limit implies, in particular, that the exponential (infinite time) 
singularity cannot be distinguished by numerical methods from the finite-time blowup with 
very large scaling coefficients aj. 

We showed that numerical results obtained by Hou and Li [261 EZ] and Grafke et al. [28] 
support the conjecture of exponential scaling of flow singularity [2^125] . The analysis shows 
that the singularity may be described by a universal self-similar solution given by the RG 
theory. The universality provides an effective criterion to be considered in future numerical 
studies. One should take into account, however, that universality may be sensitive to system 
symmetries. 
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